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The production of one hard jet in association with missing transverse energy is a major LHC 
search channel motivated by many scenarios for physics beyond the Standard Model. In scenarios 
with a weakly interacting dark matter candidate, like supersymmetry, it arises from the associated 
production of a quark partner with the dark matter agent. We present the next-to-leading order 
cross section calculation as the first application of the fully automized MadGolem package. We 
find moderate corrections to the production rate with a strongly reduced theory uncertainty. 



I. INTRODUCTION 



Since the LHC started running at a center-of-mass energy of 7 TeV searches for new physics are a major effort, 
realized in a rapidly increasing number of publications [T]. Inclusive searches for supersymmetry at the LHC have 
started to constrain the allowed parameter space of the minimal supersymmctric Standard Model [2], most notably 
in the part of the squark-gluino mass plane which can be described in terms of gravity mediation. Such searches are 
based on jet production from squark and gluino decays and two stable lightest supersymmetric particles (LSP). The 
latter could be a dark matter agent with a weak-scale mass. 

The main production mode for jets and dark matter particles at the LHC would most likely be squark or gluino pair 
production, mediated by the strong interaction [5] . The limitation of this channel is that it will be hard to extract any 
model parameters beyond the masses of the new particles [I] . The production is governed by the strong interaction 
and the (sum of) branching ratio(s) leading to jets plus missing transverse energy can be expected to be close to 
unity. Therefore, it is worth studying additional production processes which directly involve the weakly interacting 
sector of the new physics model. In supersymmetry, those are the associated production of a gluino [S] or a squark 
with a neutralino or chargino [B] 

PP -> QX° 9X ± ■ (!) 

The leading order Feynman diagrams for this process we show in Fig. [T] This channel naturally leads to one single 
hard decay jet and missing energy. This signature is not unique to supersymmetry or other models with quark partners 
and a weakly interacting dark matter agent; it also constitutes the theoretically most reliable signature for large extra 
dimensions 7\. In that sense, observing single jet production with missing energy would be one of the most exciting 
anomalies to interpret at the LHC. 

Aside from the quark-gluon and squark-gluon QCD vertices, the leading-order process is driven by the q-q-X 
interaction. Because the dominant light-flavor quarks only have a tiny Yukawa coupling, this interaction relies on 
the two weak gauge charges of the quark-squark pair involved. This way, it carries information on the composition 
of the dark matter candidate Xi an d an accurate measurement would also allow improved predictions for the direct 




Figure 1: Feynman diagrams for the associated production of a squark and a gaugino to leading order. 
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detection and relic density of dark matter. Finally, this process probes the supersymmetric relations between the 
gauge couplings and the gaugino-quark-squark couplings. Prospects for carrying out such a measurement out of the 
mono-jet rates at the LHC have been studied in Ref. |15j . 

Next-to-leading order QCD contributions to this production process arise at order aa 2 s and originate from both 
QCD (gluon mediated) and SUSY-QCD (gluino mediated) contributions. They involve virtual one-loop corrections 
as well as real emission off the initial state partons or the final state squark. While we expect the size of the virtual 
corrections combined with real emission to stay moderate, there are potential sources of huge corrections; for example, 
the production of a squark pair where one of the two squarks decays to a quark and a neutralino contributes at the 
same order ac? s , but can effectively be an overwhelming two particle production process. To avoid double counting 
and instead allow for a clear separation and prediction of the two channels, the treatment of on-shcll singularities 
in associated squark and neutralino production is crucial. Note that this separation without double counting serves 
a theoretical purpose. For an actual observable, we always need to combine QCD pair production with associated 
production because initial state jet radiation at the LHC can be as hard as decay jets [12] and cannot be distinguished 
event by event. 

Finally, we use this production process to show how heavy particle pairs at next-to-leading order can be described by 
the fully automized tool MadGolem. It generates all tree-level diagrams and the corresponding helicity amplitudes 
in the MADGRAPH framework [53], based on Helas [23]. The virtual corrections are generated with Qgraf [55] 
and numerically computed by Golem |26] , Supersymmetric counter terms are part of the model implementation. 
On-shell subtraction terms in the Prospino scheme [HI UEj can also be generated automatically. Further details 
we provide in the appendix. To our knowledge, this process is the first fully automized beyond-the-Standard-Model 
next-to-leading order (NLO) computation based on a (soon-to-be) public add-on to a major Monte Carlo generator, 
meant to be used independently by the LHC experimental community. 

II. S-UP PRODUCTION IN ASSOCIATION WITH THE LSP 

In this first section we focus on total rates for urXi and ulXi production and features characterizing the NLO 
effects. The results we present in terms of the consistent ratio K — o- NLO /cr LO , an approach which we will modify 
once we study distributions in Sec. |III| 

To avoid the current LHC bounds on squark and gluino production [5] we use a modified SPSla |20j point 
(SPSlaiooo) where the gluino mass is increased to 1 TeV. Because the gluino does not appear in the LO Feyn- 
man diagrams, its increased mass will merely reduce the impact of SUSY-QCD corrections which are nevertheless 
fully included — some effective theory issues with the proper decoupling of the gluino we discuss in detail later. 

In our numerical analysis we use the CTEQ6L1 and CTEQ6M parton densities with five flavors [TT]. For the 
strong coupling we consistently rely on the corresponding a s (fiR). Its value we compute using two-loop running 
from Aqcd to the required renormalization scale, again with five active flavors. For the central renormalization and 
factorization scales we use the average final state mass = = (m§ + m^)/2, which has been shown to lead to 
stable perturbative results [TT] [T5] . 
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Table I: Individual production rates a(pp — > qxi) and corresponding K factors for the modified SPSlaiooo scenario. The 
first and second generation squark masses happen to be degenerate. The scales are set to their central values fx R — [1% = 
(uiq + m^o)/2. In the last line we show the sum of all contributions. 
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Before we discuss the features of the NLO corrections in detail for the (dominant) ux± production channel, in Tab.[l] 
we quote the individual cross sections for all light-flavor qx\ channels. 

The results in Tab. [I] clearly reflect the flavor-locked nature of the process; ordered by the squark flavor in the final 
state, all contributions stemming from sea quarks are essentially irrelevant. Our single jet signature is driven by the 
u and d contributions, with the second generation contributing at the 5% level and, as we will see, within the NLO 
scale uncertainty. Bottom-induced sbottom production we expect to be further suppressed even though moderately 
large collinear logarithms might enhance such a signature for very light sbottoms. 

In addition, the modified SPSla parameter point is a fairly generic scenario for the weak sector and accommodates 
a relatively light mostly bino LSP. As a consequence, the neutralino coupling strength to the right and left squarks is 
substantially different, i.e. g U u LX °/g U URxl — 0-176 ~ 1/6. This explains the difference between the LO cross sections 
for ulXi an d urXi production of roughly one order of magnitude. As a bottom line, we see that for a bino LSP more 
than 80% of the total squark-LSP production rates comes from the urXi contribution. 



A. Real and virtual corrections 



In Fig. [2]we show the total cross sections as a function of the squark mass. To assess the relative impact of the real 
emission versus virtual corrections, and to spell out the differences between the right and left squarks, we only show 
results for urXi and ulXi production. The former drives the bulk of the overall qx.i event rate. The two up squark 
masses we vary simultaneously with a fixed splitting m^ L — m^ R — 20 GeV. As we can see, the virtual corrections 
are the dominant NLO effects, leading to a NLO correction of the order K ~ 1.4. 



As real corrections we consider all contributions to the NLO cross section with a three-particle final state. Virtual 
corrections include gluon and gluino induced loops, but also integrated dipoles following the Catani-Scymour dipole 
prescription |27j . While the dipole subtraction always covers the soft and collinearly divergent phase space regions, 
in terms of a variable parameter a they can be defined to extend more (a = 1) or less (a <C 1) into the non-divergent 
phase space regime |28j . Unlike for the distributions shown in Sec. Ill in this section we use a = 1, as in the original 
implementation. 

Given this choice, in Fig. [2] we see that real corrections are generally small compared to their virtual counterparts. 
They exhibit an interesting feature, namely a positive (negative) contribution for the ulXi (urXi) channel. We 
can understand this feature through the different couplings to the neutralino: real corrections to the ulXi channel 
mainly depend on g U u L ^o, but at next-to- leading order also g u ^ RX a contributes. For example, this happens for real 
corrections triggered by the gg, uu and uu initial states, as shown in Fig. [3j This way, the NLO production rates no 
longer factorize with the LO quark-squark-neutralino coupling. Such higher order effects are sensitive to the coupling 




Figure 2: Cross sections a(pp —¥ Ur/lXi) (t°P panels) and K factor (bottom panels) as a function of rriu R , L assuming 
mn L — m,u R = 20 GeV. For negative contributions to the total rate we show \a\. The remaining MSSM parameters are fixed 
to our benchmark point. Real and the virtual corrections are separated using the original Catani-Seymour dipoles [27] with 
a = 1, the integrated dipoles are included in the virtual corrections. The LHC energy is 7 TeV. 
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Figure 3: Example Feynman diagrams for real emission corrections to q\RX° production via intermediate on-shell states. 



imbalance at our parameter point, recalling g U u R xl/9uu L x1 — 6 owing to the bino-like nature of the LSP dark matter 
candidate. This is why the real corrections to the urXi channel are much smaller, and of the opposite sign than for 
the ulXi channel. 



B. Renormalization 



Before we describe the MadGolem combination of the virtual and real corrections we need to ensure that ultraviolet 
and infrared divergences are properly separated. A crucial ingredient to our higher order calculation therefore is the 
automatic treatment of the ultraviolet divergences, i.e. renormalization. In the partonic LO cross section there appear 
three parameters which require renormalization in QCD or SUSY-QCD. First, the squark mass is renormalized in 
the on-shell scheme, which means that all contributions from gluon or gluino exchange are absorbed into the mass 
counter term. For QCD corrections this scheme can be fully extended to mixing squarks [9] [21] , but in this work we 
assume the squark mixing angle to lead to negligible effects. Second, we renormalize the strong coupling constant in 
the five-flavor MS scheme and explicitly decouple the heavier colored particles from its running. This zero-momentum 
subtraction scheme |18l 1 1 9 j leaves the renormalization group running of a s to be determined only by the lightest 
colored particles. It corresponds to the definition of the measured value of the strong coupling, for example in a 
combined fit with the parton densities. The renormalization constant then reads 
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in terms of A e = 1/e — + log(47r). The coefficient of the QCD beta function we can decompose into contributions 
from light and heavy particles, (3q = /3q + /3q , with 
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The number of active flavors is n/ = 5. Correspondingly, the gluon field renormalization constant can be written as 
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Reflecting the underlying Slavnov- Taylor identities, the finite parts of both renormalization constants are related as 

S Zq = —2 S Zg g . 

Finally, we need to compensate for dimensional regularization and the MS scheme breaking supersymmetry through 
the mismatch of two gaugino and D — 2 = 2 — 2e gauge vector degrees of freedom [TBI [30] . This is done by a finite 
SUSY-restoring counter term in the quark-squark-gaugino Yukawa coupling g q q X with respect to the associated gauge 
coupling. For the q-q~X coupling, this prescription translates into a shift in = e/s w [T7]. The SUSY-restoring 
counter terms can be computed using dimensional reduction, i.e. the DR scheme. 

The fact that the q-q-x couplings has a supersymmetric limit means that we cannot numerically decouple the 
gluino contribution from this process and obtain something like a scalar leptoquark limit. While this limit is perfectly 
renormalizable, a very heavy gluino breaks supersymmetry and the ultraviolet structure of the coupling. We have 
checked that to decouple the gluinos we have to absorb a logarithm log m-g into an explicit decoupling [21] , in analogy 
toEq.[2] 
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Figure 4: NLO contributions from intermediate on-shell particles contributing to urXi production in the uu initial state. We 
show the behavior as a function of the width-to-mass ratio Fa/ma of the on-shell squark, which acts as a mathematical cutoff 
in the Prospino subtraction scheme [8], see the text for details. Fhe masses are different from our usual benchmark point, to 
illustrate all different channels. Fhe FHC energy is f4 FeV. Virtual corrections are not included. 



C. On-shell subtraction 



The Feynman diagrams in Fig. [3] naively lead to very large corrections to the Born process. The reason is that for 
on-shell intermediate states we would compare a 2 — > 2 leading order process proportional to aa s with a next-to- 
leading 2 — > 2 contribution proportional to a 2 , multiplied by a branching ratio which can be close to unity. However, 
the same diagrams are also counted for example towards qq* production, including an on-shell decay q —> qxi- To 
avoid double counting and to not artificially ruin the convergence of the perturbative QCD description of these 
production channels we remove all on-shell particle contributions from the associated production, while the same 
Feynman diagrams with off-shell propagators count towards qxi production. To provide a reliable rate prediction, we 
need to define a subtraction scheme which only removes the squared on-shell amplitudes and which does this point 
by point over the entire phase space. This Prospino scheme [8] is defined as a replacement of the Breit-Wigner 
propagator 



\m\ 2 m ^ |x| 2 Kx) |xi 2 K-) 



0(s - 4m?) Q(iriq - m^) . (5) 



In this form and for external squarks Tg is only a mathematical cutoff parameter. It can be chosen as their physical 
width (as done in the Mc@Nlo implementation 22! ) or as the uniquely defined and gauge invariant small-width 
limit 

1 7T 

lim ^—^ = — — S(s Q y - m 2 s ) . (6) 

In MadGolem this subtraction is automized. In Fig. [4] we show the numerical dependence and the stability of our 
implementation of two simultaneous on-shell subtractions, namely of uru* r and uru* l , production, where example 
Feynman diagrams are shown in Fig. [3] The intermediate gluino only appears in the uu process. For values of 
r/m < 10~ 2 the result is numerically stable, leading to the limit shown in Eq.Q. It is obvious that including the 
integrable interference terms with the continuum production is necessary to obtain a stable and numerically correct 
result. This correct description of the on-shell divergences also limits the numerical impact of higher-order corrections 
to ulXi production proportional to the 'wrong' VtUr-Xi coupling. 

Once the on-shell diagrams are properly subtracted, we can analyze the relative numerical impact of virtual and 
real corrections to UR Xi an( i Xi production in Fig. pi The different contributions are defined through dimensional 
regularization and (integrated) Catani-Seymour subtraction terms with a = 1 [27] . 
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Figure 5: Relative size of the real and virtual corrections as a function of ma R/L with a constant mass splitting ma L — m^ R 
20 GeV. Contributions from quark and squark self-energies lie below 1% and are not explicitly shown. 



Genuine QCD effects constitute the dominant contribution. The bulk of the virtual corrections arises from gluon- 
mediated qqg vertex corrections and the integrated dipoles. Each of them accounts for a +20% shift in the cross section 
for our benchmark point. The strongly suppressed gluino-mediated SUSY-QCD effects are explained by the heavy 
mass suppression, effectively decoupling the gluino from the theory. Correspondingly, the large vertex corrections 
are essentially flat as a function of the squark mass. The relative size of the integrated dipoles slightly increase with 
growing rriq towards the gluino mass range. 

Corrections to the q-q~Xi vertex can reach up to 5% for light squark masses. The box diagram, in turn, gives a 
negative contribution at the same level. Even milder is the effect from the loop-induced g-g-Xi form factor, but with 
a different sign for the ulXi 01 the urXi process. This again reflects the fact that they are sensitive to both the 
9uu L x° an d 9uu R x° couplings and hence the loop-induced g-g-Xi form factor no longer factorizes with the Born term. 



D. Scale dependences 

One of the main reasons to base LHC analyses on higher-order rate predictions is the stabilization of the dependence 
on the unphysical renormalization and the factorization scales. In Fig. [6] we present the cross sections as a function 
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Figure 6: Profile of the renormalization and factorization scale dependence for pp — > UrXi- The plot traces the scale dependence 
following a contour in the fiR-fiF plane covering fi — (0.1 — 10)/i° as shown in the left panel. We assume our benchmark parameter 
choice and VS = 7 TeV. 
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Figure 7: Total cross section for pp — > UrXi including the scale uncertainty. The band corresponds to a variation fx°/2 < 
fi R , F < 2^°. We assume our benchmark parameter choice for all masses except for ma R and show results for \/S = 7 TeV and 
14 TeV. 



of the renormalization and the factorization scales varied independently. We show the cross section profile, both at 
leading order and at next-to-leading order, moving along a contour in the hr-[1f plane. The contour is defined in 
the left panel of Fig. [6j Individually changing one of the two scales we expect a monotonous behavior, while in the 
diagonal a distinct maximum appears for some processes. Again, we only show the dominant urxI channel for the 
modified SPSla benchmark point with a heavy gluino. 

The stabilization becomes apparent in a considerable smoothing of the cr(fi) slope in the entire fiR-fJ.F plane. This 
graphical representation also indicates that the simultaneous variation of both scales would have lead to a similar 
uncertainty estimate, but this feature is clearly process and parameter dependent. Associated production cross 
sections have the feature that they are proportional to one power of a s at leading order, which means that unlike 
Drell-Yan-type channels the renormalization scale dependence is visible at leading order. In contrast to QCD pair 
production the renormalization scale does not completely dominate the combined scale dependence, though. The 
NLO uncertainty band from the scale variation ranges around (Aer)/er < 20%, down from up to 70% at leading order. 

Finally, in Fig(7]we show LO and NLO cross section predictions including the scale dependence. We see that the 
two bands are close to overlapping for fi°/2 < /j.b,,f < 2^°, i.e. our NLO result is within the LO error estimate; on 
the other hand, we also see that a conservative scale variation should not be chosen any smaller. Comparing the two 
LHC energies we see that for 14 TeV the same number of signal events corresponds to an increase in the squark mass 
by at least 200 GeV. 



E. MSSM parameter space 

In Tab. [n] we survey the predictions for squark-LSP production for all SPS points. We define an inclusive q\i 
cross section by including the first-generation squarks, q — ul, ur, d^, d,R, but without any approximation about the 
individual squark masses. As we have shown, the contributions from the second generation can be safely neglected. 
While the set of SPS benchmark points is by no means a thorough coverage of the MSSM parameter space, this 
brief scan is a useful starting point and will help us understand effects from the MSSM parameter space. The whole 
idea of the implementation of NLO corrections in the soon-to-be public MadGolem package is of course to allow for 
parameter scans by any user. 

First, we see that the K factors are largely insensitive to the specific SPS point. This follows from the dominance 
of genuine QCD effects, namely the gluon-mediated u-u-g vertex corrections and the real gluon emission. Changes in 
the squark and gluino masses do not leave a recognizable fingerprint in the relative size of the NLO corrections, which 
typically ranges around 40%. The only exception are huge corrections for ul in SPS6 and ur,cIr in SPS9, which are 
due to essentially vanishing LO rates. To next-to-leading order rate does not factorize with the LO couplings and 
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Table II: Summed cross section and corresponding K factors for all four first-generation squark processes pp — > qxl in different 
SPS benchmark scenarios. The scales are chosen at fiR p. All masses are given in GeV. 

is instead based on additional conjugate couplings. Second, we see that the total cross sections do show a strong 
correlation with the SPS points; the reason is twofold: on the one hand there is kinematics, i.e. the cross sections 
strongly depend on the final state masses in phase space; in addition, we see dynamics effects, where the strength 
of the XiQQ coupling changes substantially from one scenario to another. For relative light spectra in the range 
rriq ~ 500 — 700 GeV and m^o ~ 100 GeV, as is the case for SPSlaiooo, SPS5, and SPS6, the NLO cross sections range 

around tens (hundreds) of femtobarns at v5 = 7 (14) TeV. For TeV-scale squark masses, as in SPS2, SPS8 or SPS9, 
the cross sections stay below the femtobarn level. The K factors remain around +40%. This simple pattern hardly 
depends on the gluino mass, since the SUSY-QCD corrections are sub-leading. 

III. COMPARISON WITH MULTI-JET MERGING 

While experimental analyses based on NLO cross section incorporate significant improvements of the central values 
and the theory uncertainties, we need to ensure that this picture also includes the main distributions. From earlier 
studies we know that the transverse momentum and rapidity distributions of the heavy particles are relatively stable 
with respect to higher-order corrections [HI HB] ■ Moreover, at least for the production of heavy particles QCD jet 
radiation should be well described by the parton shower, because the collinear approximation includes sizeable prj 
relative to the masses in the final state [12]. Nevertheless, we can check quantitatively how well the NLO distributions 
from our fixed-order MadGolem computation agree with multi-jet merging jT3J [T3] - As a comparison and to estimate 
the associated theory uncertainties we use the Mlm scheme with up to two hard jets, as implemented in Madgraph. 
Any additional jets are well described by the parton shower |12) . 

An important detail when using Catani-Seymour dipoles to compute distributions is the assignment of the uninte- 
grated and integrated dipole contributions to phase space points. In the original paper, Ref. [57], it is spelled out that 
the unintegrated dipole contribution should not be counted towards their (n + l)-particle or in our case 3-particle 
phase space point, but towards the reduced 2-particle phase space. This ensures that for example the a dependence 
of the integrated and unintegrated dipoles exactly cancels not only for the total rate but also for all distributions. 
On the other hand, this implies that the pxj distribution remains unmodified from leading order and diverges at 
small transverse momenta. To turn this distribution into a useful prediction we need to consistently include a parton 
shower [36] • Obviously, neither the leading jet distributions nor the recoiling heavy system's distributions should be 
used from the fixed order computation. 
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Figure 8: Squark and neutralino pr distributions at the LHC (y/S = 14 TeV) for SPSlaiooo- We compare the merged sample 
with the fixed-order NLO computation. Both curves are normalized. We also show contributions to the NLO cross sections 
from the leading order, virtual and real parts. The latter are separated using the C at ani- Seymour dipole with a = 0.01. 



In Fig. [8] we show the transverse momenta of the squark and the neutralino as computed by MadGolem and using 
Mlm jet merging. For the Mlm simulation we cannot subtract the on-shcll singularities 1 . However, intermediate 
on-shell squarks with a subsequent decay into a neutralino will produce considerably harder neutralinos through the 
decay phase space. Therefore we only evaluate the dominant ug initial state. To the hard ug — > urXi process we add 
two hard and additional parton shower jets. 

The normalized transverse momentum distributions agree well between the next-to-leading order and the jet merging 
approaches; the slightly harder result from jet merging can be attributed to additional recoil jets. This is consistent 
with the observation that the real emission contribution to the NLO result is slightly harder than its leading order and 
virtual counterparts. In contrast to Fig. [2] where we use a = 1, these distributions are computed with a — 0.01, i.e. 
introducing the unintegrated subtraction term only very close to the soft and collinear poles. As a result, the virtual 
corrections are now smaller than the real emission contributions. Any physical observable is of course independent of 
the choice of a, as discussed above. 

The normalization of the merged Mlm sample is often an improvement over the leading order result, but not in 
a consistent quantifiable manner and to some degree dependent on the merging parameters. Therefore, it should be 
adjusted to the consistent NLO value. 



IV. SQUARK-NEUTRALINO CHANNELS 



So far, we have concentrated in the associated production of a squark with the lightest neutralino. For example 
adding leptons to the signature we would be sensitive to the production of a squark with heavier neutralinos or 
charginos. On the level of the pure production these channels often have larger rate than the bino-LSP channel. In 
Fig. M we show the total cross sections for three relatively large production channels: pp — > qj(i> qx®, QXi- The cross 
sections for both chargino charges are summed. The differences in the total cross section for fixed neutralino/chargino 
masses can be traced back to the size of the q-q~x couplings. While g U u R x° is the largest for the lightest neutralino, the 
coupling g U u L x% dominates for the next-to-lightest neutralino. The relative strength i$ U u L x%l 9uu R x\) 2 ~ 1-8 accounts 
for the bulk of the difference between the predicted cross sections. As before, we emphasize that our computation is 
based on a completely general MSSM input, encoded in a SLHA file [TO]. Negative mass eigenvalues for the neutralinos 
can be included [17]. 



1 In future releases, our PaospiNO subtraction scheme should be available for Madgraph at leading order. 
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Figure 9: Cross sections for different squark and neutralino/chargino production channels, pp — > qxi, qj(2, QXiN as a function 
of the final-state neutralino/chargino mass. We show results for \/S = 7TeV and the modified SPSla scenario. As in Tab. [n] 
we sum over all first-generation squarks. The scales are fixed to ^°n t p. 



V. SUMMARY 



We have computed the next-to-leading order QCD and SUSY-QCD corrections to the associated production of a 
squark with a gaugino. This channel could be responsible for the new physics signature of one hard (decay) jet and 
missing transverse energy. This computation is the first application of the fully automized and soon-to-be public 
MadGolem package. It makes no assumption about the supersymmetric mass spectrum, for example of light-flavor 
squarks, and relies on the SLHA interface of Madgraph. Provided its QCD counter terms are know, any new physics 
model with heavy strongly interacting particles can be included into MadGolem and its relevant processes generated 
and computed. 

For the associated production of a squark and a neutralino we find that NLO corrections due to gluon exchange 
and radiation dominate and lead to a typical correction around +40%. Gluinos only play an appreciable role when 
they are very light, which in the light of recent ATLAS and CMS measurements is becoming increasingly unlikely. 
The NLO transverse momentum distributions we have compared with a Mlm merged computation and find good 
agreement for the heavy particles produced. To avoid any dependence on merging parameters, the total cross section 
should be based on the fixed-order NLO computation. 

For associated production processes like the one considered here, it is crucial to implement a consistent and stable 
on-shell subtraction scheme separating associated production from QCD-mediated pair production. The automization 
of this subtraction in the Prospino scheme is part of MadGolem. 
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MadGolem: Automizing NLO predictions for new physics 



MadGolem completely automates the calculation of cross sections and the generation of parton-level events at 
NLO for arbitrary 2 — > 2 processes in a generic new physics framework. Its highly modular structure we illustrate in 
Fig. 10 For the generation of all tree-level Feynman diagrams and amplitudes we use Madgraph |23| . which also 
provides user interfaces and the basic code structure. The one-loop Feynman diagrams we generate with Qgraf [S3] 
and a subsequent set of specialized routines: 

1. First, we translate the Qgraf output into a code suitable for symbolic calculation languages. The structures 
describing the Feynman diagrams and the corresponding Feynman rules we rewrite as algebraic expressions, 
keeping track of external wave functions, vertex couplings and internal propagators, color factors, Lorcntz 
structure, and the overall sign from external fermion fields. With this modification MadGolem can cope with 
genuine features of new physics processes, such as Majorana fermions. 

2. Then, we map the analytical evaluation of the color, helicity and tensor structures onto partial amplitudes, i.e. 
a basis of color, helicity and tensor structures based on the spinor-helicity formalism. 

3. In addition, we apply an analytical reduction to scalar loop integrals, based on a modified Passarino-Veltman 
reduction scheme implemented in Golem [2*5] . 

4. Finally, we combine the results with the ultraviolet counter terms (which are also generated automatically). The 
final output for the virtual corrections we return both as analytical Maple and numerical Fortran90 code. 

To compute complete NLO rates MadGolem uses dedicated modules for the different types of divergences ap- 
pearing in NLO calculations. To remove infrared divergences from real and virtual gluon emission we resort to the 
Catani-Seymour dipole subtraction [27] . Our automatic dipole generation is a modified and extended version of Mad- 
Dipole now including all massive Standard Model and supersymmetric dipoles. All integrated dipoles explicitly 
retain the dependence on the phase space coverage of the subtraction term a [28] . Intermediate on-shell divergences 
we remove following the Prospino scheme described in Sec. II C 8J. Its local subtraction terms we generate automat- 



ically. The ultraviolet renormalization counter terms we generate with Madgraph amplitude. They are expressed in 
terms of two-point functions which are supplied in a separate library. The current library supports QCD corrections 
for arbitrary 2 — ► 2 processes in the Standard Model and the MSSM and can easily be generalized to other new 
physics models. 



MadGolem does not require the user to interfere with the computation of NLO cross section and parton-level 
events from the basic input. The model and the process are specified through Madgraph input cards. Options like 
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multi-particle notation are supported together with additional specifications that allow us, for instance, to separate 
QCD from SUSY-QCD contributions or retain subsets of one-loop contributions. 

We have performed exhaustive checks to ensure the reliability of MadGolem. Total NLO cross sections we have 
tested both in the SM and the MSSM and covering numerous initial/final states, interactions, and topologies. Cancel- 
lation of all divergences and gauge invariance of the overall result have been confirmed numerically and analytically. 
The finite renormalized one-loop amplitudes we have compared with FeynArts, FormCalc and LoopTools re- 
sults |32| . Particular attention we have paid to numerical stability and convergence of our Catani-Seymour dipoles 
and the on-shell subtraction. Finally, we have checked our final results with the literature (e.g. e + e~ — > q*q [33]) as 
well as with Prospino [3J] and MadFks [35] (e.g. pp -> 1*1). 
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